Resistance spikes and domain wall loops in Ising quantum Hall ferromagnets 
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We explain the recent observation of resistance spikes and hysteretic transport properties in Ising 
quantum Hall ferromagnets in terms of the unique physics of their domain walls. Self-consistent 
RPA/Hartree-Fock theory is applied to microscopically determine properties of the ground state 
and domain-wall excitations. In these systems domain wall loops support one- dimensional electron 
systems with an effective mass comparable to the bare electron mass and may carry charge. Our 
theory is able to account quantitatively for the experimental Ising critical temperature and to explain 
characteristics of the resistive hysteresis loops. 
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For two-dimensional electrons in a perpendicular mag- 
netic field B±, independent electron eigenstates occur 
in manifolds known as Landau levels with macroscopic 
degeneracy AB±/§o, where A is the sample area and 
$o is the magnetic flux quantum. The zero-width en- 
ergy bands are responsible for a tremendous variety of 
many-body physics that has been observed in the quan- 
tum Hall regime [QJ^). Quantum Hall ferromagnetism, 
of interest here, occurs when two different Landau levels 
distinguished by the cyclotron energy, spin, or quantum 
well subband labels of their orbitals are brought into en- 
ergetic alignment and the Landau level filling factor v 
is close to an integer. Neglecting charge fluctuations, 
low-energy states of quantum Hall ferromagnets (QHFs) 
are specified by assigning to each orbital in the Landau 
level a two-component spinor (cos 9/2, e t(p sin 9/2) corre- 
sponding to a pseudospin oriented along a general unit 
vector fh = (sin cos 95, sin # sine/?, cos 9). (The influence 
of remote Landau levels can be captured perturbatively 
as necessary.) While ordered states can occur when any 
two Landau levels simultaneously approach the chemical 
potential, the nature of the ground state is sensitive to 
the microscopic character of the crossing Landau levels 
|,| . Isotropic II, XY |-| and Ising QHFs §§£1 
are now well established. Our work is motivated by the 
recent observation |k| of hysteretic transport and un- 
explained resistance spikes when Landau levels with dif- 
ferent quantized kinetic (cyclotron) energies cross. We 
argue that the resistance spikes are due to charge trans- 
port in the ID quasiparticle systems of long domain wall 
loops and establish a correspondence between their oc- 
currence and vanishing domain-wall free-energy density 
at the Ising transition temperature, T c . 

The dependence of the uniform QHF state energy 
per electron on pseudospin orientation has the form [Q: 
E[fh] — —bm z — Jm 2 z , where b is an effective magnetic 
field that includes both single-particle Landau level split- 
ting and interaction contributions H and J > is an 



effective Ising interaction parameter. At b = 0, the 
m z = 1 (pseudospin f) and m z = — 1 (pseudospin |) 
states are degenerate. In the following we establish an 
association between b = and the experimental resis- 
tance spikes, and propose an explanation for the spike 
origin. For tilted magnetic fields and variable 2D electron 
densities, the b = condition at a given filling factor is 
achieved along a continuous line in the two-dimensional 
(Btot — B±) space, which can be explored experimentally 
by tilting the field away from the sample normal. (Btot 
is the total magnetic field.) 
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FIG. 1. Main graph: Theoretical (dashed) line of the 
n — 2 spin-up and n = spin-down level coincidences in the 
perpendicular magnetic field - total field space, calculated us- 
ing nominal sample parameters of De Poortere et al. [[l3j and 
taking into account RPA-screening of the Coulomb interac- 
tion. The density of the 2D electron system is 2.5 x 10 11 cm -2 . 
The solid line in the main graph is a fit to the experimental 
magnetoresistance spike positions presented in Ref. in 
Fig. 3(A). Inset: Effective external field of the Ising model 
calculated as a function of the perpendicular magnetic field 
at a fixed field tilt angle 34°. Landau level configurations at 
positive and negative values of the 6-field are also indicated 
in the inset. 
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In Fig. H we compare our theoretical |||l4| 6 = line 
for v = 3, based on numerical self-consistent-field calcu- 
lations for the geometry of De Poortere et al.'s sample 
and on many-body RPA/Hartree-Fock theory with the 
line along which resistance spikes were observed. The co- 
incidence of these two curves strongly suggests that the 
spikes occur when 6 = 0. The same calculations 0] yield 
the estimate J = 0.0l8e 2 /e£/ ks ~ 2K, where £ is the 
magnetic length defined by 2ir£ 2 B± = $o- 

The ground state of an Ising QHF has m z — 1 for 
6 > and m z = — 1 for 6 < 0. At finite temperatures, 
non-trivial pseudospin magnetization configurations be- 
come important. For Ising QHFs, an elementary calcu- 
lation shows that spin- wave collective excitations have a 
gap A/fcs = AJ/ks ~ 8K. Since the hysteretic resis- 
tance spikes occur only for T < 0.5K, spin-wave excita- 
tions cannot play a role. Instead, as we now explain, the 
important thermal fluctuations in Ising QHFs involve do- 
main walls between m z — 1 and m z — — 1 regions of the 
sample. 

In classical 2D Ising models the critical temperature 
can be understood as a competition between unfavor- 
able near-neighbor-spin interaction energy along a do- 
main wall, Irf, and the wall configurational entropy, 
Ls c = L/^fcflln(3), where £ is the domain wall persis- 
tence length. Both give free-energy contributions propor- 
tional to wall length, L, with the former effect favoring 
short walls and the latter contribution, which is propor- 
tional to temperature, favoring long walls. For T > T c 
the system free energy is lowered when domain walls 
expand to the sample perimeters, destroying magnetic 
order. The structure of domain walls is more compli- 
cated in Ising QHFs. As the domain wall is transversed, 
the local pseudospin orientation goes from the north pole 
(m z = 1) to the south pole (m z = —1), at a fixed orien- 
tation ip of its x — y plane projection. We have evaluated 
the energy per unit length 7 of an infinite domain wall 
by solving self-consistent Hartree-Fock equations: 



tanflpO 



H^(X) 



H ],] (A) 



(i) 



where the Hartree energy is given by 

H^,{X) = J2j d 3 n J cPriVin -f 2 ) x 
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and the exchange energy by 
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In Eqs. (H) and (||), V{?\ — r 2 ) is the RPA-screened 
Coulomb potential and the self-consistent-field one- 
particle orbitals, ip a> x{r)-, are extended along the domain 
wall and localized near wavevector k dependent guiding 
centers X = k£ 2 . The energy density 7 is proportional to 
the increase in Hartree-Fock quasiparticle energies, inte- 
grated across the domain wall. We find that the domain 
wall width is typically several magnetic lengths and for 
the v = 3 coincidence we find that j£ = 0.009e 2 /e£. 

A unique property of QHFs is the proportionality be- 
tween electron charge density and pseudospin topologi- 
cal index density. It is this property that is responsible 
for the fascinating skyrmion physics extensively studied 
in the isotropic case JIH. In the case of Ising QHFs, 
the proportionality implies a local excess charge per unit 
length along a domain wall pu — eVtp ■ n/(27r) where 
n specifies the local direction along the domain wall. 
Single- valuedness of the magnetization requires that the 
winding number of the angle if around a domain wall 
loop be quantized in units of 2ir and hence that the ex- 
cess charge of a domain wall loop be quantized in units 
of the the electron charge e. The free-energy associated 
with the classical tp field fluctuations within a domain 
wall, 



U = jk B T\nZ: Z = J Vpexp (E c [p]/k B T) 



(4) 



is controlled by the Coulomb interaction energy E c due 
to the consequent charge fluctuations. 

Assuming a domain wall persistence length £ w £, the 
free energy density of Ising QHF domain wall loops is 
given by / = /cs7Tn(3)/£ + 7 + f v and equals zero at 
T = T c . For the v = 3 QHF, these considerations im- 
ply that infinitely long domain walls proliferate and or- 
der is lost for T larger than the transition temperature 
T c w 500 mK. The close correspondence between this 
T c estimate, and the maximum temperature (430 mK) 
at which hysteretic resistance spikes are observed |pj 
strongly supports our contention that the unusual trans- 
port phenomena are a consequence of the existence of 
long domain wall loops in these materials. In the follow- 
ing we first discuss the 6-field, Landau level filling factor, 
and temperature dependence of the system's domain-wall 
soup and then demonstrate that this picture can account 
for many details of the transport observations. 

Domain wall loops are characterized by their length 
and by their charge, with infinitely long loops appearing 
only for 6 = and T > T c , For finite size loops with 
a typical radius larger than the domain wall width we 
can use our Hartree-Fock self-consistent results for 9(X) 
to estimate the Coulomb self-interaction energy. A two- 
dimensional charge density of a circular loop with excess 
charge e distributed uniformly along the domain wall is 
given by 
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Since the corresponding Coulomb self-interaction energy 
is proportional to the square of the charge and approxi- 
mately inversely proportional to the length, charged do- 
main wall loops have a higher energy and, at integer 
filling factors, will always be less common than neutral 
domain wall loops. Resistance spikes are generically ob- 
served slightly away from integer filling factors, however, 
and here the situation changes because the domain wall 
loops can exchange charge with the rest of the 2D elec- 
tron system. 
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FIG. 2. Near-Tc, v > 3, statistical weights of domain wall 
loops with excess charge from — e to 2e are plotted as a func- 
tion of the loop radius, a) For bjks ~ 20mK, the effective 
external field corresponds to the edge of the measured resis- 
tance spike. Small neutral domain wall loops dominate in this 
case, b) For b « 0, the coincidence condition is approached 
corresponding to the spike maximum. The typical domain 
wall loops are large and carry an excess charge. 

The lowest energy elementary charged excitations of 
the v = 3 ground state are ordinary Hartree-Fock elec- 
tron and hole quasiparticles jl7| , not charged domain wall 
loops. In systems with no disorder, the chemical poten- 
tial lies in the middle of the Hartree-Fock gap when v 
is an integer but moves quickly (by 8fi) toward the elec- 
tron quasiparticle energy for v > 3 and toward the hole 
quasiparticle energy for v < 3. These chemical potential 
shifts will be reduced by disorder which broadens the 



quasiparticle bands. The change in chemical potential 
favors charge Q over neutral skyrmions by a large factor 
exp(Q|5/x||/fcsT). This factor can be estimated quanti- 
tatively using the experimental value of the quasiparticle 
excitation gap (~ 2 K) 

Also important in controlling the domain wall soup is 
the effective field b, which measures the distance from 
Landau level coincidence. For b ^ the energy of a 
domain wall loop has a contribution proportional to b 
and the number of condensate electrons contained within 
the loop: 

E b = ~^2 J drr [cos 9(r)-l] . (6) 

This contribution will decrease the number of large do- 
main wall loops enclosing the minority phase and is inde- 
pendent of the charge carried by the loop. Summing up 
E},, the Coulomb self-interaction and chemical potential 
contributions and the Hartree-Fock domain wall energy 
of a circular loop of radius R, 2irR , y, we can estimate sta- 
tistical weights of neutral and charged domain wall loops 
in the sample of De Poortere et al. In Fig. [2] we plot our 
results for temperature near T c and v > 3. For non-zero 
6-fields, small neutral domain wall loops dominate, while 
typical loops at the coincidence are large and carry an 
excess charge. 

We now address characteristic features of the mea- 
sured resistive hysteresis loop. Dissipation can occur in 
Ising QHFs as a result of Hartree-Fock quasiparticle dif- 
fusion, charged domain-wall-loop diffusion, or as a result 
of charge diffusion within domain- wall loops. It is clear 
that the resistance spikes, which appear only for small b 
and T < T c , are associated with the appearance in the 
sample of large domain-wall loops. Even though these 
loops tend to be charged at the spike maximum, we ex- 
pect that they will be immobile because of their large size 
and that dissipation due to their motion is small. Instead 
we propose that mobile quasiparticles inside domain walls 
are responsible for the increase of dissipation. In Fig. ^ 
we plot the Hartree-Fock quasiparticle energies in a cross- 
section of a domain wall, obtained from Eqs. In 
the center of the domain wall the quasiparticle gap is re- 
duced by nearly 50%. Away from integer filling factors, 
the bottom of these ID quasiparticle bands will lie below 
the chemical potential which is pinned to the bulk quasi- 
particle energies. We note that, unlike quantum Hall 
edge states, counter propagating states exist within each 
loop. At v s» 3, for example, the ID states have a nearly 
parabolic dispersion characterized by an effective mass 
to* « 2m e . The particles can cross the sample by scat- 
tering between overlapping loops. For v — 3 and T w T c , 
it follows from inset of Fig. |l| and from Fig. |^ that the 
characteristic loop radius at the spike edge is 3£. At low 
temperatures or small magnetic lengths (high 2D elec- 
tron gas densities), domain wall loops become small and 
dilute, loops do not overlap, and charge diffusion within 
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domain walls cannot contribute to dissipation. This ex- 
plains the absence of resistance spikes at Landau level 
coincidence under these circumstances. 




FIG. 3. Hartree-Fock quasiparticle energies, measured 
from the integer filling factor chemical potential, are plot- 
ted as a function of orbit guiding center for zero effective 
external field b. For \X\ ^> the spectrum describes ordi- 
nary, bulk electron and hole quasiparticles. Near the domain 
wall center a parabolicly dispersing ID channel is formed. 
The thin and thick lines are obtained assuming pseudospin-up 
and pseudospin-down majority polarizations, respectively. At 
v = 3 (solid lines), the screening of the Coulomb interaction 
is substantially stronger in the latter case resulting in dif- 
ferent quasiparticle spectra for the two configurations. The 
inset illustrates schematically the corresponding peak asym- 
metry in the measured magnetoresistance near filling factor 
v = 3. For comparison, we also show in the main graph spec- 
tra for v = 5 (dashed lines). The negligible difference between 
pseudospin-up and pseudospin-down majority polarizations is 
consistent with the similar peak heights recorded in up and 
down magnetic field sweeps near v — 5. 

The above mechanism also explains the different resis- 
tance spike heights observed in up and down held sweeps 
fbif near v = 3. As shown in the inset of Fig. [j], the ma- 
jority pseudospin Landau level for the up-sweep is the 
n = 2 spin-up level while for the down-sweep it is the 
n = spin-down Landau level. This difference in Lan- 
dau level configurations in the two sweep directions alters 
remote Landau level screening in the sample which has a 
marked effect on the quasiparticle energy spectrum. The 
reduction of the quasiparticle gap in the domain wall, rel- 
ative to its bulk value, is stronger in the up-sweep case, 
leading to more domain- wall quasiparticles and more dis- 
sipation, as seen in experiment. Similar agreement be- 
tween the measured hysteresis loop properties and do- 
main wall quasiparticle spectra applies for v — 4 |l6| ] . In 
Fig. H we also plot energy spectra at v = 5 which are 
nearly identical for the up or down majority pseudospin 
orientations. This explains the absence of peak-height 



asymmetry in the hysteresis measurement at this filling 
factor v = 5. |ll| 
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